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The observation of B — > rp, + /if at present experiments would be a clear sign of new physics. 
In this paper we calculate this process in a 2HDM framework where the decay is mediated by the 
exchange of spin zero particle with flavour changing neutral current couplings. If we identify the 
scalar with the the newly discovered state at LHC with a mass ~ 125 GeV then we get BR(B a — > 
rfl + fif) ~ 10~ 6 and BR(Bd —}Tfl + fif) ~ 1CP 7 . We also calculate this process in minimal 
supersymmetric standard model and find BR(B S — > rfi + fif) ~ 10 -11 . 
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I. INTRODUCTION 



Flavor changing neutral current (FCNC) events are rare in the Standard Model (SM) both in the quark and 
the lepton sectors. These processes can be affected by new physics (NP). FCNC involving the third family 
quark and leptons are interesting as their larger masses make them more susceptible to new physics effects in 
various extensions of the SM. As an example, in certain versions of the two Higgs doublet models (2HDM) 
the couplings of the new Higgs bosons are proportional to the masses and so new physics effects are more 
pronounced for the heavier generations. Moreover, the constraints on new physics involving, specially the third 
generation leptons and quarks, are somewhat weaker allowing for larger new physics effects. Hence, in the down 
quark sector there is intense interest in the search for NP in FCNC b — > s and b — > d transitions in Bd, s mixing, 
b — > s(d)qq and b — > s(d)£ + £ ~ decays. Many measurements in the B sector have constrained NP effects though 
slight discrepancies from SM predictions still remain in several decays !lj . In the up quark sector one can look 
for NP is the processes t — > c(u) transitions Q. 

In the lepton sector FCNC decays are severely suppressed and lepton number for each family is conserved 
in the SM. Therefore, FCNC effects and lepton flavor violation in this sector are sensitive probes of NP. The 
presence of neutrino masses and mixing suggest NP in the lepton sector which in turn could enhance FCNC 
effects and lepton flavor violation. It is interesting that measurements involving the the tau lepton show 
some discrepancies from SM expectations. The branching ratio of B — > tv t shows some tension with the SM 
predictions 3 though more recent measurement are consistent with the SM Q • This could indicate NP 0] , 
possibly coming from an extended scalar or gauge sector. There is also a seeming violation of universality in 
the tau lepton coupling to the W suggested by the LEP II data which could indicate new physics associated 
with the third generation lepton Q. 

More recently, the BaBar collaboration with their full data sample of an integrated luminosity 426 fb" 1 has 
reported the measurements of the quantities [t| 

= BR{B_^Dt- V - t ) = 
K ' BR(B -> Dl-vt) 

R(D*) = BR ^ B - ~* D * T ~ Vr ] = 0.332 ± 0.024 ± 0.018 . (1) 
v ' BR(B ->■ B*l-v t ) y ' 

The SM predictions for R(D) and R(D*) are 0-Q 

R(D) = 0.297 ±0.017, 

R(D*) = 0.252 ±0.003, (2) 

which deviate from the BaBar measurements by 2a and 2.7er respectively. The BaBar collaboration themselves 
reported a 3.4cr deviation from SM when the two measurements of Eq. [1] are taken together. This result, if 
confirmed, could indicate NP involving new charged scalar or vector boson states [l3|. Models with new charged 
scalars or charged vector bosons also include neutral particles that can cause FCNC effects. 

In this paper we will focus on the process B s d — > iilj . This involves two FCNC interactions- one in the 
quark and one in the lepton sector. Typically one believes that this process will be very small even with new 
physics and in many NP models this is indeed true. This is because in these NP models this FCNC process 
arises in loops and is suppressed. In this paper we will focus on the decay B — > r/2 + (if and make predictions 
for its rate. This decay was considered in earlier papers in various extensions of the SM [ll[. In this work we 
will calculate the process in two scenarios representing two extension of the SM. 

In the first scenario we will consider this process mediated by a scalar which we will call X in a two higgs 
doublet framework (2HDM). We will identify X with the newly discovered particle, with mass ~ 125 GeV, 
observed at LHC [12|,ll3| with supporting evidence for its existence from Fermilab 14|. In this case B — > Tp+(if 
is generated at tree level due to the presence of bqX and t/iX couplings. The coupling bqX can be constrained 
by Bd,s mixing while the T(iX couplings can be constrained from rare r decays. Recently in Ref. [l5, 16| these 
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coupling were considered and it was noticed that the t/iX coupling could be of similar size as the SM higgs to 
tt coupling. Since B s ^ are pseudoscalars then clearly X must be a pseudoscalar or have both parity conserving 
and parity violating couplings to quarks and leptons. The spin parity of the X particle is not yet known though 
various strategies that allow for the spin and parity determination are being actively pursued [17]. In this 
scenario the branching ratio for B — >■ rp + /if can be significant and can potentially be observed at LHCb. In 
the second scenarios we will consider a popular extension of the SM- the minimal supersymmetric SM (MSSM). 
Here the FCNC couplings arise via loops and so our expectation is that the branching ratio for B — > rp + /if 
will be tiny. The main point here is that if this process were to be observed at the LHCb in the near future this 
would indicate the presence of a rather non-traditional NP with appreciable tree level FCNC couplings in the 
lepton sector and a light spin zero mediator. 

The paper is organized in the following manner. In the next section we describe the general effective Hamil- 
tonian for b — > sl^lj decays. In the following section we consider the decay B — > rp + /if in a 2HDM framework 
with FCNC higgs couplings. This is followed by a calculation of B — > rp + /if in MSSM. Finally we summarize 
our results and present our conclusions. 



II. EFFECTIVE HAMILTONIAN 



In this section we discuss the effective Hamiltonian for b —> q{= d, s)^ lj transitions. We will focus on the 
b — > sl~lj Hamiltonian as the b — > dl~lj~ Hamiltonian can be obtained from it with the obvious replacements. 
The effective Hamiltonian for the quark-level transition b — » sl^lj (I = e, /i, r) in the SM is given by 



H 



SM 
eS 



4G F 

' V2 



V, 



ZVtb {Y, <?i(/-0<W + C 7 [sa^(m s P L + m b P R )b] 



i=l 



c 9 { - srPL b) L% + c 10 2s (s^P L b) }, 



(3) 



where = k^f^lj, — ^7^75^-, and Plm — (1 T 7s)/2- The operators Oi (i — 1, ..6) correspond to the Pi 



in 18j, and m b = mb{n) is the running 6-quark mass in the MS scheme. We use the SM Wilson coefficients as 
given in fl9j |. 

The total effective Hamiltonian for b — » sl~l^~ in the presence of new physics operators with all the possible 
Lorentz structure excluding tensor once can be expressed as [2(3] 

n cS (b -> alrlf) = H^i 4 + U V J + U S J , (4) 

where V.^ 1 is given by Eq. ([3]), and the NP contributions are 

U V J = -N F [r v (srfP L b) L% + R A (s 7 »P L b) 1% + R' v {s^P R b) L% + R' A (*y»P R b) L%} , (5) 
U S J = —Np {r s (sP R b) L i:i + R P (sP R b) L% + R' s (sP L b) L V] + R' P {sP L b) i|} , (6) 

where Np = V t * V t b, Lij = Ulj, and = k^lj. In the above expressions, Ri, R'^i = V, A, S, P) are 

the NP effective couplings which are in general complex. 

In terms of all the Wilson coefficients, the branching ratio for the decay B s — > Ijlj can be obtained from 



Br(B s ^lri r) = ^^ mB ^ V ^ 



* I 2 

ts\ 



64tt 3 



rrii + rrij ' 
m Ba 

2 



{ I 1 - P^^r 1 ) 1 |Fy(m * - mj) + Fs|2 + 1 1 - (^^r 1 ) ] lFA{m * + mj) + Fp|2 }' (7) 



m B s 
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where 

Fy = Cg + Ry — Ry 7 

Fa = Cio + Ra — R'a > 
F s = r x (R s ~ R' s ) , 

F P = r x (R P - R' P ) , (8) 

m 2 

with r x — n^^+ra ( M ) ■ The Wilson coefficient Fy does not contribute in the lepton flavor conserving decays 
B s — > and B s — » t~t + . Within the SM, Fa = C\q only contributes to B s — > and B s — > t~t + . 

The lepton flavor violating decay B s — > T[i is not allowed in the SM. Thus, the Wilslon coefficients Cg and Cio 
do not contribute to this decay. Next, we look at two models of NP that contribute to B s ,d — > l^lj decays. 



III. FLAVOR VIOLATING HIGGS DECAYS 



In this section we consider the model where a neutral scalar boson X with mass ~ 125 GeV and flavor- 
violating couplings [15|, [l6| mediates the B s ^ — > decays . In this model, the effective Lagrangian which 



describes the possible flavor- violating couplings of X to SM fermion pairs in the mass basis is, 15j, |l6[ : 

Cy = -mJlP R - Yy (fif R )X + h.c+ , (9) 

where ellipses denote nonrenormalizable couplings involving more than on Higgs field operator, fi, = 
are SU(2) L doublets, fp — updp,vp,lp are the weak singlets. The indices run over generations and fermion 
flavors with summation implicitly understood. In the SM the Higgs coupling are diagonal, Yij = (rrii/v)5ij with 
v = 246 GeV, but in general NP models the structure of the Yij can be different. The couplings in Eq. © will 
also lead to the decay b — > sl~lf at tree level, mediated by a virtual X . The lagrangian gives 



R 



S(P) 



1 Y sb {Y^ T +Y n 



4m x 



N F 

1 Y b * s (Y^-Y^) 



R 



V(A) 



R[ 



V(A) 



Am 2 x 
0. 



The branching ratio for the B s 



~t + decay can be obtained from Eq. ([7]) as, 



(10) 



Br(B s h~t^ 



T Bs 



fB s m B s 



h 5127r(mfc(/i) + m s (^)) 2 



mp 3 



mp. 



\Y sb n. 



Also, the branching ratio for the B s —> = /i, r) decays is, 



Br(B s rr) 



t Bs fBG F m Bs a 2 em \Vtb V t * s \ 
h 64tt 3 




4m? 



V m a ) 



Re[Y u ](Y sb -Y* 



2m x NF(mb{[i) + m s (/i)) 



Vtor / 



(12) 



where we assume Yu is real and set it to its SM value Yu — mi/ v. 

Flavor violating Higgs coupling in Eq. ([9]) can generate flavor changing neutral currents at tree level. The 
couplings Y s b and Yi, s are constrained by the mass difference AM S in the B s — B s mixing. The AB = 2 weak 
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Hamiltonian for this process can be found in the Ref. [l^ . The theoretical expression for AM, is given in 
Ref.Q. It is found that one can reproduce the measured value of AM S — (17.719 ± 0.043)ps 1 |22j if 



\Y xh - Y, 



ba 



10" 



(13) 



The leptoninc couplings Y" Mr and Y TfJ , are constrained by the decays r — > \v) and r — ¥ 3/i, and the magnetic (Sa^) 
and electric dipole (ei M ) moments. The theoretical expressions for T(t — > ^7), T(r — > 3u), Aa» = a^f p — a^ M , 
and cL are taken from The experimental bounds for these observables are taken from [23l - l25f : 



Br(r -> fij) < 4.4 x 10" 8 
Br(r 3/i) < 2.1 x 1CT 8 

Aa M = (2.87 ± 0.63 ± 0.49)10~ 9 
-10 x lQ- 20 e cm < d M < 8 x 10~ 20 e cm. 



(14) 



It is found that one can satisfy the above four experimental bounds if \Y^ T \ < 0.064 , \Y TfJ ,\ < 0.061 for mx = 125 
GeV. The fit results allow us to estimate the decay rate Br(B s — > fir) = Br(B s — > h~t + ) + Br(B s — > t~^ + ). 

The variations of the decay rates Br(B s — ► ^t) with \Y s b — F b * | and {Y^r + Y* \ are shown in Fig. 1. Our 
result predict Br(B s ^ can be as large as ~ 6 x 10" 6 for m x = 125 GeV and f B . = 0.229 ± 0.006 Q. 
In the calculation we have used the experimental constraint \Y^ T + Y*\ < 0.10. In this model, we obtain the 
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FIG. 1. The variation of Br(B s — > ^i~t + + /i + t~) with the couplings \Y s j, — Y^ 3 \ and |Y MT — Y*^\ for mx = 125 GeV 
and /b s = .229 ± 0.006 GeV. Scatter points are allowed by AMj; xp and experimental bounds in Eq. [1$ . 



branching ratio Br(B s — > \x /i + ) < 4.3 x 10 9 using Y^ — m^/v with v = 246 GeV, mx — 125 GeV and 



3.17 x 10~ 9 . Our results for Br(B s 



= (3.2tH) x 10~ 9 at 95% C.L. in Ref. 



271. Also, 



m T /v and mx — 125 GeV. The values 



f Bs = 0.229 ± 0.006 GeV. The corresponding SM prediction is 
is consistent with the current upper limit Br(B s 

we obtain the branching ratio Br(B s — > t~t + ) < 8.1 x 10~ 7 using Y TT 
of this branching ratio in the SM is « 6.81 x 10~ 7 . The latest LHCb-measurement of Td/T s implies a limit of 
Br(B s — » t~t + ) < 3% In passing we note that, the large width difference in the B s meson system can 
changes these results by ~ O(10%) as pointed out in 3oT- 32| . 



Finally, in this model we obtain the tau longitudinal polarization fraction for the decays B s 



Pl = 



Br 1 



- T+ [\ T = 1/2] — Br 1 



- T+ [A r = -1/2] _ [7 _ , 

- r+ [A T = -i/2] yL v 



m B . 



m Bs 



Br B ^- T +[X T = l/2] + Sr Bs _ 

( 2 (|y Air | 2 - iv r ;i 2 ))/{ [1 - {^^-y\ iv + n;i 2 + [1 - i 1111 ^)'} iv - n;i 2 } ■ (15) 
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The result indicates Pr. does not depend on the b — > s couplings. Fig. [5] shows the dependency of Pl on the 
quantity (|V | 2 - IF^I 2 ). 

The couplings Ydb and Ym for the decay Br(B — > /it) = BriB — > ^~t + ) + Br(B — > t~/j + ) are similarly 
constrained by the mass difference AM,j = (0.507 ± 0.004)ps _1 [23]. The measured AM^ can be reproduced if 
\Ydb — Y b * d \ ~ 10~ 4 . Using this result, we found the decay rate Br(B — > jj,r) < 2 x 10~ 7 . 



IV. SUPERSYMMETRY 



The supersymmetry contributions to the process B s — > /i + r~ are given by the one-loop box diagrams, shown 
in Fig [3l where charginos and neutralinos are exchanged. The interaction of charginos, neutralinos can be 
written as 
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FIG. 3. One loop SUSY contributions to B s -> t~/U 
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he 
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(16) 
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where the mixing matrices in the super-CKM basis are given by 



H Ul - a 
u Uij — 9 



-V; i {T Ul )j 1 + V* 2 {T u «). h 



V2M W sin /3 



V 



ji > 



m d i 



\f2mw cos (3 
mi 

\/2mw cos f3 ' 



/ 



7 D L 
A Iik 



9 

V2 
9 

V2 
9 

V2 

g r2 



(N* k2 + tan^iV* 1 )(r i -) /l - N* 3 (T L «) n 



mi 



mw cos /3 



2t<m6 w N kl (r L x) Ii +N k3 (T Ll <) Ii 



mi 



mw cos j3 



1 



m d 



V2 



(N* 2 + - tan 6 w N* kl ) (T D - ) u + N£ 3 (T D ») H 

i M\y COS p 

m d 



t- A n6 w N kl {T D «) H + N k3 (r D t) Ll 
3 mw cos [3 



VjT m , (17) 
(18) 
(19) 
(20) 
(21) 
(22) 



where /, J = 1..3 and i = 1..6. The matrices T's are the matrices that diagonalize the squark and slepton mass 
matrices, while U and V are the unitary matrices that diagonalize the charged mass matrices. Finally, N is the 
matrix that diagonalizes the neutralino mass matrix. 

One can easily show that the lightest stop gives the dominant effect for the chargino contribution. In this 
case, the amplitude is given by 



GFCCem,m 3 B /b 

A^± = — -2 V ts V tb S{m Ii ,m^,m~+,rn~+) 



/ 1 (m T + to„) 2 - F s 
m 2 R mi, + m s 



'1 



K - m^) 2 F'p - F P 
m 2 B mi, + m s 
(23) 



where 
F s 



F's 



sin 9w cos 2 /3 



mi 



V2V kl u^(((^ UL *)3 3 ) 2 + ((r^*)36) 2 ) 
nJh-.ur' <■ !,: i (T' '-') + (r % *) 36 (r^*)26)], 

(V2v kl u k2 ((T u -*) 33 (r u -) 33 + (r^*) 36 (r^)36) 



mw sin (3 
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sin 2 9w cos 2 (3 



mi 



mw sin (3 



v k2 u k2 ((r u -) 33 (r u «*) 23 + (r u -) 36 (r u «*)26) 



(25) 
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where i,j,k,l = 1..2 denote the heavy and light sparticles (stop-quark, sneutrino and chargino). 
In deriving these expression, the following hadronic matrix elements have been assumed: 



(0 |s 7 5 6| B s ) = -if Bs 



' m b + m s 

Similarly, the neutralino contributions can be written as 



(0\sb\B s ) = 0. 



(29) 



GFOt em m% Jb 
A x° = ^ 5 S ( m l i > TO ? 3 > m xl > m x1 ) 



(ra T + m^) 2 R' s — R 



m b + Tn 




(m r - m^) 2 R' p - R P 



m b + m s 



where 



D fy*D R yD L i y*D R yD L \( y*L R yL L , y*L L yL R \ 

n S — \^23k Zj 33k + Zj 2Qk ^ 36k) (^331 ^231^^331 ^231 ) 

ryl I y*D L yD R | y*D L yD R \( y*L L yh R . y*L R yL L \ 

n S ~ \ Zj 23k ^33k + ^26k Zj 36k)\ Zj 33l ^231 + ^331 ^231 ) 

j) (y*D R yD L 

Up — {^23k ^33k 



y*D R yD L \ I y*L L yL R 

■°26k Zj 36k)\ Zj 33l ^231 
7 *D L yD R 



y*L R yLL \ 

^331 ^231 ) 



p' ( y*D L yD R _|_ y*D L yD R \/ y*L L yL R 

n P — \ Zj 23k n 33k "T z '26fe Zj 36k)\ Zj 33l ^231 



y*L R yL L 
^331 ^231 



(30) 



(31) 
(32) 
(33) 
(34) 



However, it turns out that the neutralino contribution is typically smaller than the chargino effect. Also in case 
of diagonal slcpton mass matrix, which is severely constrained by the lepton flavor violation process such as 
fi — >■ e7 and r — > /ry, the neutralino contribution identically vanishes, it means that Rs = R's = Rp = R'p — 0. 
In this case, the branching ratio of B s — > h + t~ is given by 

|2 



TBa Gjai m m%J Bs \V*V tb \< 
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in 



B, 



64tt 3 
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(l-(^i 



m Bs 



-) 2 )(l-(- 



m B 



ml 



F's-Fs 



m b + m s 



(35) 



For lightest stop and lightest sneutrino masses of order 200 GeV and lightest chargino mass of order 120 GeV, 
one finds that 



S{m^.,mn m~+,m~+) ~ 1.23 x 10" 



(36) 



Also with the largest possible off-diagonal elements in the stop and sneutrino mass matrices, one may obtain: 

F s - F s > ^ 93.35, F P - F P , ~ -82.87. (37) 

Therefore, one concludes that the SUSY contribution to the process B s — > u + t~ is of order 

Br(B s -> h + t-) ~ 3.43 x 10^ n , (38) 

which is about 5 order of magnitudes smaller that the tree level contribution of two-Higgs doublet model, 
discussed in the previous section. 

V. SUMMARY 



In this paper we analyzed the rare decay B — > rp, + ui? in two possible extensions of the SM. In particular, we 
considered the tree level contribution in 2HDM, where the decay is mediated by a light scalar which is identified 
by the recently observed state at LHC. In the second example, the one loop contribution to this process in 
MSSM was computed. We found that in the former scenario the branching ratio B s — > rp + [if is of order 10~ 6 
( Bd — > Tp + [if <~ 10~ 7 ), while in the later scenario BR(B S — > /i + r~) ~ 10~ n . Therefore, one concludes that 
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probing this lepton flavor violating process at the LHCb would be a significant hint for a non-traditional new 
physics beyond the SM that has a tree level FCNC. 
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